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The transport of energetic electrons is sensitive to magnetic perturbations. By using 3D 
numerical simulation of test particle drift orbits we show that the transport of untrapped 
electrons through an open region with magnetic perturbations cannot be described by a 
diffusive process. Based on our test particle simulations, we propose a model that leads 
to an exponential loss of particles. 
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1. Introduction 

There is an increasing interest in the effect of magnetic perturbations on transport in 
fusion plasmas, as it has been shown that externally generated magnetic perturbations 
influence many aspects, including e.g. the presence and properties of Edge Localized 
Modes (Evans et al. 2004; Suttrop et al. 2011), plasma rotation (Fietz et al. 2015), Neo¬ 
classical Tearing Modes (Koslowski et al. 2006) or losses of relativistic electrons (Lehnen 
et al. 2008). Electrons with relativistic velocities (runaway electrons) are the particles 
that are most sensitive to magnetic perturbations as they follow open field lines with 
v || c. The purpose of this work is to provide insights into the effects of magnetic field 
perturbations on the transport and loss of energetic electrons by using three-dimensional 
(3D) numerical modelling of the electron drift orbits. 

Runaway electrons are frequently generated in tokamak disruptions and can cause 
damage on plasma facing components. One of the envisaged mitigation strategies is to 
use magnetic perturbations to lower the confinement of the runaway electrons and pre¬ 
vent the formation of a potentially harmful runaway electron beam. Resonant magnetic 
perturbations were tested with some success on various tokamaks (Lehnen et al. 2008; 
Yoshino et al. 1999; Commaux et al. 2011). Possible effects and uses of 3D magnetic field 
perturbations in ITER have been explored (Shimada et al. 2007; Papp et al. 2011 a,6, 
2012). The individual particle orbits in perturbed magnetic fields are chaotic (Papp et al. 
2012), and for a precise calculation of the transport one needs to resort to 3D numerical 
simulations at large computational expense. However, for the purposes of self-consistent 
modelling of runaway electron dynamics, it is important to describe the transport due to 
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magnetic perturbations with a reasonably simple model, so that it can be included in the 
modelling of the coupled dynamics of the runaway electron current and resistive diffusion 
of the electric field (Papp et al. 2013). This work aims to provide insight into how such a 
simple model can be constructed so that it faithfully reproduces the dynamics of particle 
losses due to magnetic perturbations. 

2. Effect of magnetic perturbations 

The transport of electrons in enclosed stochastic magnetic fields is usually described 
with the Rechester-Rosenbluth diffusion coefficient Drr = Tiqv\\R (SB/B) 2 (Rechester 
& Rosenbluth 1978). Here q is the safety factor, v\\ ~ c is the parallel velocity, R is 
the major radius and SB/B = (( SB/B) 2 )^ is the flux surface averaged normalized 
magnetic perturbation amplitude as a function of radius. In mixed magnetic topologies 
consisting of magnetic islands, intact toroidal magnetic surfaces and stochastic regions, 
this diffusion coefficient is not valid (Myra & Catto 1992; Mynick & Strachan 1981). 
In such circumstances, some fraction of the electrons will be trapped in magnetic is¬ 
lands and the confinement time of these particles can equal the characteristic system 
evolution time. Furthermore, in open chaotic systems, such as the edge region of a mag¬ 
netically confined plasma in the presence of external magnetic perturbations, not only the 
Rechester-Rosenbluth diffusion is not applicable, but the diffusive approach itself breaks 
down. Unlike in the case for enclosed stochastic fields, the region of increased transport 
(compared to the unperturbed case) extends all the way to the wall, making it easier for 
particles to get lost. As was confirmed by numerical simulations described in this paper, 
the runaway electron flux starting from a certain radial position is not proportional to the 
runaway density gradient, but the local runaway density itself. If a flat runaway profile 
is initialized, losses occur immediately all around the simulation zone, i.e. the flux is not 
gradient-driven, as the inital density gradient is zero. This dynamics is due to the open 
ergodic region that is created by the RMP coils. Note, that the diffusive approach may 
still be applicable in closed ergodic regions. The net radial displacement of the particles 
is highly sensitive to the initial position and can differ by several orders of magnitude. 
Particles born within resonant loss regions will reach the outer edge of the plasma faster 
than the rest and be lost to the wall. Since the system is open at the plasma edge, there is 
no return flux of fast electrons, and hence the density gradient is not the governing factor 
in the evolution, as it would be in a diffusive picture. As we will show, this process can 
be approximated with a position dependent exponential loss of particles, that describes 
the ensemble behaviour with reasonable accuracy. 

In an axisymmetric tokamak equilibrium the magnetic topology consists of nested 
magnetic surfaces that provide good particle confinement. The main loss mechanism for 
fast electrons is due to the drift orbit shift associated with the energy gain (Papp et al. 
2011 a,6), which is determined by the electric field evolution. In the simulations presented 
here we use static externally-imposed magnetic fluctuations on a static equilibrium field 
with nested toroidal magnetic surfaces. Whenever a magnetic perturbation is present, 
the topology of the magnetic field is strongly modified. There appear island chains (the 
locations of which are correlated with the low order rational values of the rotational 
transform -b = 1 /g), together with intact, but deformed magnetic surfaces (KAM barriers 
Kolmogorov (1957); Arnold (1963); Moser (1962)) at irrational -b values. In between these 
ergodic zones are formed. 

We solve the relativistic gyro-averaged equations of motion for test energetic electrons 
in the perturbed field with the ANTS code (Papp et al 2011a). The simulations have 
been carried out for a specific ITER-like scenario as described in (Papp et al. 20116). The 
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Figure 1. Electron loss time as function of initial coordinates in normalized flux (ip) and 
poloidal angle ($) at 60 kA coil current. 


perturbations are externally imposed by the 9x3 quasi-rectangular perturbation coils at 
the outer side of the device. In the configuration we use (formerly labeled “B” (Papp et al. 
2011 6)), electrical currents flowing in these coils generate a 120° rotationally symmetric 
perturbation field that decays radially towards the inside of the plasma and results in 
the relative perturbation strength SB/B > 0(1O -3 ) outside the ip = 0.5 flux surface 
(the normalized radius relates to the normalized flux as r/a ~ y/ip, where a is the minor 
radius of the torus). 

To determine the evolution of the particle ensemble in such a topology a direct calcu¬ 
lation of test particle trajectories has been carried out. In each simulation a minimum of 
40 000 test particles with 1 MeV energy were launched parallel to the magnetic field lines 
with various starting position distributions. In the simulations we use a time-dependent 
electric field obtained for an ITER-like disruption scenario calculated with a model of the 
coupled dynamics of the evolution of the radial profile of the current density (including 
the runaways) and the resistive diffusion of the electric field (Smith et al. 2009). The 
electric field peaks between 6-15 ms, which means that most of the acceleration happens 
in this time frame, as it is shown in Fig. la of Ref. (Papp et al. 2012). The electrons are 
gaining energy due to the electric field mainly during this time-interval. In the unper¬ 
turbed scenario, the drift orbit shift associated with the energy gain dictates the particle 
losses, which in this ITER case happens on the r ^ 10 ms timescale (Papp et al. 20116). 

Figure 1 shows the loss time of electrons launched with (ip G [0.5,0.8], $ e[ 0,27r], (p= 0) 
as a function of initial coordinates in the (p — 0 poloidal cross section. The simulation 
was performed for 100 ms. The particles corresponding to the black dots have not been 
lost during this time, as these were trapped in the remnant O points of island chains 
formed around low-order rational surfaces. Figure 1 clearly reveals the chaotic nature of 
the process and highlights the extreme sensitivity of the individual particle loss times to 
the initial position. 

Despite this complexity, the ensemble behaviour is smooth. Figure 2a shows the frac¬ 
tion of lost particles as a function of time for particles launched at individual flux surfaces. 
All these curves show an exponential decay of the confined particles up to the satura¬ 
tion value, determined by the fraction of particles trapped in remnant O-points of an 
island chain intersecting the starting surface. Dashed lines show exponential fits and the 
characteristic loss times are indicated above the curves. 
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Figure 2. (a) Fraction of lost particles for different starting radii as a function of time. Dashed 
lines show Ni ost (t) = Nq(1 — exp{— t/r}) fits, (b) Fraction of lost particles within ^ E [0.5,0.8] 
for different perturbation currents as a function of time. 

The magnetic perturbation enhanced net radial transport leads to particle losses on 
the sub-millisecond timescale. Our energy scans confirmed that on this short timescale 
the behaviour is similar regardless of the energy of the particles if E > 1 MeV, as the 
particles mainly follow the open magnetic field lines with the speed v\\ ~ c. The long-term 
(> 10 ms) behaviour is determined by the remnant island structures. Particles starting 
in specific spatial positions can be trapped in the remnant O-points and confined even 
at energies of 100 MeV. 

If the magnetic perturbation strength is increased (by increasing the current in the 
coils) the island structure will change non-linearly. Therefore the relationship between 
the time-evolution of the losses and the magnetic perturbation magnitude is not trivial. 
Figure 2b illustrates the fraction of lost particles as a function of time for three different 
values of the perturbation current. To highlight the effect of the magnetic perturbation 
magnitude, in these simulations the electric field has been set to zero (the electrons are 
not accelerated as they were in the other simulations). 

The smooth ensemble behaviour is also reflected in the evolution of the electron den¬ 
sity as a function of radial position (t/>) and time, as is shown in Fig. 3a. The continuous 
radial loss picture is only distorted by the enhanced confinement of the remnant island 
chains at ^ =0.5, 0.6 and 0.7. For comparison, Fig. 3b shows the evolution of density due 
to the Rechester-Rosenbluth diffusion for the same initial conditions. The evolution of 
the density n(^, t) is largely different: the diffusive model predicts better confinement at 
the edge but worse confinement in the core, by orders of magnitudes. The generation of 
runaway electrons is exponentially sensitive to the existing amount of runaways through 
the avalanche process. If instead of a test particle scenario the runaway electron gener¬ 
ation is included in the calculation, the differences between the two cases are increased 
even further, as will be shown later in this paper. 

Following the observations illustrated in figures 2a-b and 3a, we approximated the 
particle losses with an N\ ost (t) oc 1 — exp{—trend, where r(^) is the characteristic 
loss time associated with a certain initial *0 position for the particles. We have performed 
fits of r(0) to the evolution of fast electron density obtained by the 3D test particle 
simulations for various cases, such as the ones illustrated in figures 2a-b and 3a. The 
fitted values of r(0) are shown in figure 4b for three different perturbation magnitudes. 
We note that the dependence of r(0) on the radial coordinate is relatively smooth and 
is well approximated by an exponential dependence r(0) oc exp{— 'ip/'ipo} (dashed lines), 
where 0o oc 5B/B. Local maxima in r(0) correspond to remnant islands and KAM zones 
in the perturbed magnetic topology. 
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Figure 3. (a) Evolution of fast electron density as calculated by the ANTS code, 

(b) Evolution of the density due to Rechester-Rosenbluth diffusion. Both figures are calculated 
for 60 kA coil current. 



Figure 4. (a) Evolution of electron density using the fitted r( 'ip) values for the 60 kA configura¬ 
tion, including particles confined in the remnant islands ( n c ). (b) Fitted values of r as a function 
of normalized flux for 3 different perturbation amplitudes. 


Replacing the diffusive transport in a numerical calculation with the exponential loss 
model, we can calculate the evolution of fast particle density using the fitted r(val¬ 
ues. A pure exponential loss description cannot reproduce the confined fraction. How¬ 
ever, if we split the particle density into a confined and non-confined fraction n(^, t) = 
n c (V0 + ^ncC*/M); where the confined part is constant at a given radial position (ex¬ 
tracted from the ANTS simulations) and the non-confined part follows the exponential 
loss n nc (£, ip) oc 1 — exp{—£/r(t/>)}, we can reproduce the results of the 3D simulation. 
Figure 4a shows a calculation based on the exponential losses for the 60 kA case, includ¬ 
ing the n c ( , 0) particles confined in the remnant islands. It clearly resembles the particle 
evolution calculated by ANTS (figure 3a - the figures are placed next to each other for 
easier comparison). 


3. Self-consistent calculations 

To illustrate the differences between the implications of the diffusive and exponential- 
loss models, we have implemented both loss models in a self-consistent simulation of 
the runaway current and electric field evolution (Papp et al. 2013). Figure 5 shows a 
disruption simulation for the same ITER scenario that was used for the 3D simulations. 
The thermal quench time is set to 0.2 ms and the effective charge Z e ^ = 2 to mimic the 
parameter evolution of massive gas injection simulations (Hoffmann et al. 2015; Papp 
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Figure 5. Effect of diffusive (solid lines) vs exponential (dashed lines) losses on plasma (RE) 
current evolution for ITER in a self-consistent calculation for different perturbation strengths. 

et al. 2013). As a general conclusion for all cases, following the thermal collapse, at 
^3 ms virtually all of the Ohmic current is replaced by runaways, and the runaway 
current consists of ~10% seed and ^90% avalanche current. Without losses a runaway 
beam of ^9 MA would form. Using the Rechester-Rosenbluth diffusion model and the 
5B/B profile caused by the perturbation coils we get suppression of the runaways after 
20 ms at /rmp = 60 kA and after 10 ms if /rmp = 120 kA. /rmp = 30 kA not only does 
not lead to suppression, but through the redistribution of the primary seed current leads 
to temporary stronger avalanching, increasing the overall runaway current. However, as 
we argue in this paper, with the support of the 3D test particle calculations, the losses in 
these systems are not well described by diffusion. Also, the simulations show that although 
the perturbation penetrates to the core, the transport is essentially untouched inside the 
radius r/a = 0.7 for /rmp = 60 kA. Therefore, even if the transport is largely increased 
towards the edge region in the exponential loss description, this does not contribute much 
to runaway suppression. The reason is that most of the runaways are formed in the plasma 
core, where there is good runaway confinement according to the test particle simulations. 
As the dashed lines in figure 5 clearly illustrate, virtually no runaway suppression is 
achieved even at /rmp = 60 kA with the exponential loss model. At /rmp = 120 kA we 
see a slow drop of the runaway current, however, /rmp > 75 kA will not be accessible in 
ITER due to hardware limitations. 


4. Conclusions 

The most important conclusion of the present work is that the transport of the un¬ 
trapped electrons in open chaotic fields cannot be described by a diffusive model. Instead 
there is an exponential decay of confined particles, where the characteristic loss time 
is an exponential function of radius and the perturbation strength. This difference has 
a decisive effect on the runaway current and the associated electric field evolution in a 
self-consistent calculation. While diffusive losses would suggest that runaway avalanches 
can be avoided by using strong enough externally applied magnetic perturbations, the 
exponential loss model shows that this does not have to be the case. Although it is 
difficult to predict exactly the magnitude of the loss-time and its radial dependence in 
a specific scenario without 3D numerical simulations, the exponential dependence and 
non-diffusive nature of the transport is robust. Therefore implementation of such a model 
in self-consistent fluid simulations should give important insights into runaway electron 
dynamics and the prospects of RE mitigation in the presence of magnetic perturbations. 
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